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ABSTRACT 

We present MAMBO 1.2 mm observations of five BzA'-pre-selected vigorous starburst galaxies at z ~ 2. 

Two of these were detected at more than 99.5% confidence levels, with 1.2 mm fluxes around 1.5 mJy. These 
millimeter fluxes imply vigorous activity with star-formation rates (SERs) «5OO-15OOM0 yr^\ confirmed 
also by detections at 24/rm with the MIPS camera on board of the Spitzer satellite. The two detected galaxies 
are the ones in the sample with the highest SERs estimated from the rest-frame UV, and their far-IR- and UV- 
derived SERs agree reasonably well. This is different from local ULIRGs and high-z submm/mm selected 
galaxies for which the UV is reported to underestimate SERs by factors of 10-100, but similar to the average 
BzX-ULIRG galaxy at z ~ 2. The two galaxies detected at 1.2 mm are brighter in K than the typical NIR- 
counterparts of MAMBO and SCUBA sources, implying also a significantly different X-band to submm/mm 
flux ratio. This suggests a scenario in which z ~ 2 galaxies, after their rapid (sub)mm brightest phase opaque to 
optical/UV light, evolve into a longer lasting phase of X-band bright and massive objects. Targeting the most 
UV active BzKs could yield substantial detection rates at submm/mm wavelengths. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — galaxies: starburst 
— cosmology: observations — submillimeter 


1. INTRODUCTION 

In the last few years, selection techniques at different wave¬ 
lengths have revealed a variety of, at first glance, different 
high redshift source populations (e.g., Steidel et al. 1996; 
2004; Smail et al. 1997; Kurk et al. 2000; Eranx et al. 2003; 
Daddi et al. 2004a; Cimatti et al. 2004). Highly debated ques¬ 
tions are how these objects are linked to each other; how much 
they overlap with each other; and what the evolutionary path 
is, if any, between them and with the different types of low 
redshift galax ies. Observations in the local univers e are in¬ 
dicating te.g.. lThomas et al~ll2005tlNelan et alJl2005h that the 
stars in massive early-type galaxies formed in general at high 
redshift, z >> 1, with short formation timescales, suggesting 
high peak star-formation rates (SERs). Populations of high 
redshift galaxies, among those mentioned above, having large 
stellar masses and high SERs have indeed also been found. 
However a detailed characterization of the formation of early- 
type galaxi es is still mis s ing. 

Recentlv. lDaddi et alJil2004bll presented a technique based 
on optical/near-IR photometry in the B, z and K bands that, 
at least to Kvega < 20, allows to obtain virtually complete 
samples of galaxies at redshift 1.4 < z < 2.5, including star- 
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forming and passively evolving galaxies. This technique is 
reddening independent for z ^2 star-forming galaxies, that 
in X < 20 samples have SERs estimated between 50 — 1000 
solar masses per year. These SERs require large and pos¬ 
sibly uncertain reddening corrections. The average long- 
wavelength emission properties of BzK star-forming galax¬ 
ies (sBz Ks hereafter) at z = 2 have been measured accu¬ 
rately bv iDaddi et al.^^2005^ ■ finding that the typical sBzK is 
an Ultra Luminous Infra-Red Galaxy (ULIRG) with SFR ~ 
200 — 300 M© As the BzK technique is substantially 

complete for z ~ 2 galaxies, one expect that also submm/mm 
selected galaxies (SMGs, see Blain et al. 2002 for a review) 
shoul d fulfill the BzK crite ria, as many of them lie at z 2 
(e.g.. IChanman et^l2005h . and that the most active sBzK 
galaxies could actually be SMGs. In this letter we present 
observations with the 117- element Max-Planc k Millimeter 
Bolometer Array (MAMBO: iKrevsa et al.ll9^ of five lumi¬ 
nous {Kvega 19) sBzK galaxies, with high expected SERs, 
based on their optical/UV properties. We report the detection 
of two of the five galaxies, and discuss some implications on 
the relation between sBzK and SMGs, and on the assembly 
process of massive galaxies at high redshifts. We assume a 
Salpeter initial mass function (IMF) from 0.1 to 100 M©, and 
a cosmology with — 0.73.0.27, and h = //o[km s^* 

Mpc-*]/100 = 0.71 ifSnergel et ali2003L 

2. SAMPLE SELECTION 


We selected for MAMBO observations five bright Kvega < 
19.2 mag yBzK-galaxies (see Table [D , taken from the 
700 arcmin^ wid e ’Daddi-Field’ at 14h (iDaddi et alJl200(i 
iKong et ^120051) . All the targets had spectroscopic redshift 
identifications from deep VLT-VIMOS optical spectroscopy 
and sometimes VLT-SINFONI and/or Subaru-OHS/CISCO 
near-IR spectroscopy, as a part of our o ngoing survey for 
BzK selected galaxies dKong et alJ 120051) . The optical and 
near-IR spectroscopic observations will be presented else¬ 
where. Redshifts range from 1.6 to 2.8, and SERs estimated 
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TABLE 1 

Properties of the sBzK galaxies observed with MAMBO. 


Source 

(1) 

RA 

(J2000) 

(2) 

DEC 

(J2000) 

(3) 

K 

(AB) 

(4) 

BzK 

(AB) 

(5) 

Zspec 

(6) 

E(B-y) 

(7) 

SFRyv 

(M@lyr) 

(8) 

Mass 
(10*'Me) 

(9) 

Si .2 mm 
(mJy) 

(10) 

SFR„ 

(MQ/yr) 

(11) 

S 24 ^m 
(Ally) 
(12) 

Comments 

(13) 

OBJ4415 

14:48:35.35 

08:52:27.0 

20.92 

1.61 

2.53 

0.49 

177/140 

2.8 

0.55±0.54 



Lya; AGN-2 

OBJ4193“ 

14:48:33.23 

08:54:14.1 

20.67 

1.99 

2.76 

0.40 

~200 

4.0 

-0.26±0.58 



Lyc(; AGN-2 

OBJ2742 

14:49:20.50 

08:50:52.3 

20.84 

0.02 

2.16'’ 

0.53 

666/400 

1.4 

1.34±0.43^ 

700-1400 

70±23 

UV-abs/IR-em? 

OBJ2426 

14:49:28.96 

08:51:52.9 

20.76 

0.05 

2.37 

0.59 

546/1400 

1.8 

1.50±0.42 

800-1900 

200±50 

UV-abs/IR-em 

OBJ 1901 

14:49:41.38 

08:59:50.5 

20.58 

0.68 

1.60 

0.33 

186/140 

1.7 

0.54±0.43 


79 ±25 

UV-abs/IR-em 


Note. — Col. (1 ): Object identifica tion. Cols. (2)-(3): Optical/NIR coordinates. Col. (4): Total ^-band magnitudes. Col. (5): BzK-co\or. BzK = 
[z — K)ab — ^ z)ab IDaddi et all20n4bh . Col. (6): Spectroscopic redshift. Col. (7): Colour excess: £"(5 — E) — 0.25[(S — +0.1], Daddi et al (2004b). 

Col. (8): Two SFR estimates are from the UV corrected 1500 A luminosities IDaddi et^l2004bD and from SED fitting, respectively. Col. (9): Stellar mass 
estimated from SED fitting; uncertainties are typically a factor of two. Col. (10): MAMBO “on-off” flux at 1.2 mm. Col. (11): SFRs estimate based on L[r. 
Col. (12): MIPS 24 (im flux. Col. (13): Comments to the UV/optical spectra of individual sources. OBJ4193 is at the edge of the Subaru imaging data, 
photometry is less accurate. this redshift is lower quality, although supported by UV-spectrum and possible Ha detection. Excluding the first observing 
consisting of 7 scans under mediocre weather conditions would give S \,2 mm — 1.57 + 0.46 mJy (T 1.2 mm ^ 0-2). 


from the dust corrected UV luminosity range from ^ 150 to 
^ 1400 M© Although the K-band brightest sBzK natu¬ 
rally tend to have fairly red average optical/UV colors, prefer¬ 
ence among sources with measured redshifts was given to ob¬ 
jects with red B — z colors, corresponding to large E{B — V). 
The reasons for this preference are twofold: 1) optical red¬ 
ness might be an evidence for the presence of large amounts 
of dust, enhancing the chance of detections with MAMBO; 
2) the red colors en sure these objects wo uld rarely qualify for 
the z = 2 criteria of ISteidel et al.l ll2004h . thus allowing to in¬ 
vestigate to which extent vigorous starburst galaxies may be 
missed by surveys selected in the UV. 

The field has b een observed in the X-ray with XMM- 
Newton for 80 ks JBrusa et alJl2005h and with Chandra as a 
mosaic of 3 x 30 ks with an overlapping area of 90 ks (Brusa 
et al. in preparation). Targets for MAMBO observations were 
chosen to be un-detected in the X-rays (corresponding to X- 
ray luminosities ^ 10^^ erg s^*) and with no strong AGN 
feature in their UV/optical spectra. However, two of the tar¬ 
gets (OBJ4415 and 4193) show faint and narrow type2 AGN 
emission lines in their (almost continuu m-less) UV spectra, 
simila r to some of the objects described bv ivan Dokkiim et alJ 
ll2003l) . These two objects were not detected at 1.2 mm. 

3. OBSERVATIONS AND RESULTS 

Each of the five sBzKs were observed under predominantly 
excellent weather on several days during December 2004 with 
MAMBO-117, at the IRAM 30 m telescope on Pico Veleta. 
MAMBO operates at an effective frequency of 250 GHz, cor¬ 
responding to 1.2 mm (FWHM w 10.7 The science tar¬ 
gets were observed in “on-off” observing mode, which is the 
photometric mode of MAMBO, for 2700 to 5700 s. This ob¬ 
serving technique is based on the chop-nod technique where 
the target is placed on a reference bolometer element (on- 
target channel). A detailed description of t his now stan¬ 
dard o b serving technique can be found in e.g.. lBertoldi et al] 
( l2000l) ;lLutzet 70 il2005l) . In order to reduce systematic er¬ 
rors and avoid spurious detections, the targets were observed 
several times and on different days. The data were reduced us- 
ing M OPSIC*^. an upgraded version of MOPSI ilZvlka et alJ 
Il998h . Every scan with its individual subscans were care¬ 
fully inspected for the presence of possible outliers or anoma- 

http://www.astro.ruhr-uni-bochum.de/nielbock/simba/mopsic.pdf 


lies or influence of high opacity (Ti ,2 mm > 0.2 — 0.3). Erom 
each channel we subtracted the correlated skynoise from the 
surrounding channels. The outstanding weather conditions - 
opacity below T 1.2 mm < 0.2 during 60% of the time - ensured 
stability and high data quality over the several days in which 
the data where gathered. 

A flux of about 1.3-1.5 mJy (significant at 3.2 — 3.6(7 level) 
is detected in the reference bolometers corresponding to the 
positions of two of the five sBzK galaxies (OBJ2742 and 
OBJ2426; see Table [TJ. The MAMBO array of detectors 
consist of 117 bolometers spread over a sky-region of about 
4 arcmin in diameter. In order to establish reliably the sig¬ 
nificance of the detections, and to estimate the likelihood 
that these may be just be obtained by chance, we indepen¬ 
dently reduced the observations for each of the 117 MAMBO 
bolometers. Considering the 5 independent objects observed, 
this provides a large control sample of about 550 reliable 
bolom eters. Given the source counts at 1.2 mm (iGreve et al J 
ITonl . the chance probability to observe a galaxy with flux 
Z 1.5mJy within half a beam FWHM is 0.6%. However, we 
assume conservatively that all the control bolometers are see¬ 
ing a blank mm sky. No control bolometer has a detection 
with S/N > 3.6 (OBJ2426), while 3 control bolometers have 
S/N > 3.2 (OBJ2742). Therefore we conclude that the de¬ 
tections are reliable at better than the 99.5% confidence level. 
Fig .[Oshows beam sized fields in the B and z band around the 
two detections. 

In order to constrain better the IR-SED of both MAMBO 
detections, and to further improve the confidence in the detec¬ 
tions, we obtained MIPS 24/rm imaging data of the field from 
the Spitzer archive. In the 500s per sky-pixel data, OBJ2742 is 
detected at a 3 (7 level with a flux Sza fim ~ 70 jiJy. OBJ2426 
is brighter with a flux S24 fim ~ 200 ± 50 jxJy (the large er¬ 
ror being due to the fact that this object is blended with a 
nearby source 3.5" to the west, see Fig. Q] that has a sim¬ 
ilar 24/rm luminosity). The 24 ixmll .2 mm flux ratios for 
our two sBzKs are fully consistent with the range reported by 
llvison etTI] ^2004^ for MAMBO galaxies, and also consisten t 
with SCUBA sources dFraver et al.l2004llEgami et alJ2004h . 

The two galaxies detected are the ones with the highest es¬ 
timated SFRs based on the rest-frame UV. No other clearcut 
differences e.g., color, magnitude, could be found between 
the MAMBO detected and undetected sBzK galaxies. A pos¬ 
teriori, it is reasonable and reassuring that we detect with 
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Fig. 1.— B-band (top) and z-band (bottom) images for the fields around 
the two MAMBO detections (dashed circle): OBJ2742 and OBJ2426. The 
beam size of MAMBO (10.7" FWHM, big circle) is fully contained in the 
images. Both sources show an in'egular, possibly merger-like moiphology. 
North is up and east to the left. 


MAMBO the 2 sources that were estimated already to be the 
most active in our sample. This further enhance the conh- 
dence that the detections are reliable. 

4. DISCUSSION 

The fluxes measured at 1.2mm (close to the IR-peak of the 
dust emission) and at 24/rm with Spitzer (in the mid-IR rest- 
frame region) imply that the two detected sBzKs are active 
IR galaxies. To obtain accurate estimates of their total IR 
luminosities (8-1000/4m) would require further observations 
between 24/im and 1.2mm, and especially in the range 60- 
lOO/rm rest-frame, that could constrain the temperature of the 
dust. For z = 2, obtaining these data to the required depths is 
beyond the possibilities of the available observing capabili¬ 
ties. However, we can derive constraints on the total IR lu¬ 
minosities (and thus SFRs) by assuming these 2 detections 
have IR SEDs similar to those of well studied z = 2 popula¬ 
tions, i.e. the average sBzK galaxies (Daddi et al. 2005) or 
SCUBA selected galaxies (Chapman et al. 2005), or by as¬ 
suming that local empirical correlations between mid-IR or 
submm luminosities and Lir (e.g. Chary & Elbaz 2001) holds 
also at z = 2. We And that this suggests IR luminosities in 
the range of Lir w 3-10 x IO'^Lq , which correspond to SERs 
in the range 5OO-15OOM0 the factor of 3 uncertainty 
being indicative of the range of estimates from the different 
methods. 

We considered the possibility that the detected mm fluxes 
are due to AGN emission. The submm-to-X-ray spectral 
slope asx (Alexander et al. 2003) is an indicator of AGN 
activity. Eor the two MAMBO-detected sBzKs we derive an 
upper limit of ayy 1 .20*^ that, if compared to Eig. 7 in 
[Alexander et alJ ( 120031) . points toward the starburst part of 

We assume ratios of 2.5 — 3 between 850pm and 1.2 mm fluxes, that 
are the factors spanned by the models at the luminosities of our objects. 
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Fig. 2.— (Top Panel) SFR/^ vs. SFRf/u for the five sBzKs observed by 
MAMBO (red open squares). For compaiison, w e show the median SFR of 
sBzKs from the GOODS-N field (blue filled square;|Da ddi_etalJ2QQ3), a sam¬ 
ple o f local ULIRGs (open triangles: IGoldadere^ni20Q?Trrentham et all 
II99^ . and the median for radio-identified SMGs from the SC UBA redshift 
survey (open diamondi lSmail et alJ2Q04tlChaDman et all2QQ5ll . The dashed 
lines represent different ratios of SFR/^ to SFRf/v from 1 to 100. (Bottom 
Panel) K-band magnitude (Vega) vs. 1.2 mm fiux for our two sBzKs de¬ 
tected by MAMBO (red open squares). For comparison, we show the me¬ 
dian of sBzKs from the GOODS-N field (blue filled squa re), radio-identified 
MAM BO galaxies from the NTT Deep Field (crosses; IDannerbauer et all 
120041) . and SMGs with ^ > 1.5 from the SCUBA redshift survey (open di¬ 
amonds). The dashed lines represent different ratios of 1.2 mm to K-band 
flux from 10 to 2000. 


' 2000 X 

looo"'". + 

+ 


200 'A o 


o 


100 

50 


A>|, 


o UAo 

- O + -xy 






40 


+ - 


■ - 10 


1 .2mm flux [mJyl 


the diagram occupied, e.g., by SCUBA galaxies in HDFN 
and Arp 220, and exclude both luminous unobscured (3C273) 
and obscured quasars. Of course, this does not allow us 
to completely discard the presence of Compton thick AGNs 
(NGC 6240). However, the lack of any AGN features in the 
UV spectra, together with the reasonable agreement between 
UV- and far-IR- estimated SERs, suggest that the far-IR activ¬ 
ity of these sources is powered by star formation. 

An interesting result of our observations is indeed that the 
IR luminosities and related SERs inferred for our objects from 
the mm fluxes are in quite a good agreement, at most a factor 
2-3 larger than those estimated from the dust-reddening cor¬ 
rected UV luminos ities (Fig.|2l top pa nel). This agreement, 
already noticed by iDaddi et'Sl ll200-5t) for the typical sBzK 
galaxy, is somewhat surprising. A fairly different situation 
seemingly holds for submm/mm selected galaxies (Chapman 
et al. 2005; Swinbank et al. 2004) and for local far-IR se¬ 
lected ULIRGs (Goldader et al. 2002, with only a notable 
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exception), where the observed IR luminosities and far-IR es¬ 
timated SFRs exceed by 1-2 orders of magnitude the same 
quantities estimated from the dust corrected UV luminosities 
(Fig. 121 top panel). Extending this comparison further (Fig.|2l 
bottom panel), we notice that sBzKs galaxies in this paper and 
in Daddi et al. (2005) are generally brighter in the K band (~ 
18 — 20 mag) than the typical counterparts of MAMBO g alax- 
ies ( iBertoldi et alj|20^ |Dannerbauer et al.||2002[ |2004]) and 
SCUBA sources (e.g.. iBlain et ali2002l:ISmail et alJ2004) . al¬ 
though they lie at a similar redshift z ~ 2. The mm/submm 
fluxes of SMGs are instead generally higher than those of 
sBzKs (Fig. Ell. This implies a significantly different /T-band 
to submm/mm flux ratio between far-IR selected vs. optically 
selected z ^ 2 galaxies. 

Similarities also exist between sBzKs and far-IR selected 
populations both in their irregular morphologies (Conselice et 
al. 2004; Trentham et al. 1999; Daddi et al. 2004a; this paper, 
see Fig.^ and in their red observed UV slopes (i.e. high ap¬ 
parent reddening; Chapman et al. 2005; Daddi et al. 2004b; 
2005; Goldader et al. 2002; this paper, see Table 1). The dif¬ 
ferences (submm/mm flux levels, fir-band to submm/mm ra¬ 
tios and SFRuv to SFRjr ratios) and similarities (morpholo¬ 
gies, UV slopes) are most likely effects induced by the differ¬ 
ent selection of the samples (fT-band versus far-IR), and can 
be reasonably understood by assuming that the far-IR emitting 
region of most SMGs is observed behind a thick dust screen, 
nearly opaque to the optical and near-IR radiation. This can be 
fit into an evolutionary scenario in which the starburst galax¬ 
ies are initially deeply embedded into a large column density 
of dust, during the submm/mm brightest phase when substan¬ 
tial obscuration of a few magnitudes or more is present even 
at observed //-band (rest-frame 5000 — 8000 A for z = 2 — 3). 
The UV light is almost completely attenuated behind the dust 
screen so that it is impossible to measure the ongoing SFR 
from the residual UV light transmitted. As time passes by. 
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SFR and Ljr are reduced and the dust is diffused and ejected 
out of the main star-formation regions (presumably by galac¬ 
tic winds or other kinds of feedback), or simply destroyed. 
Thus the resulting galaxy is brighter in the near-IR and trans¬ 
parent enough in the rest frame UV that the transmitted light 
can account for the ongoing SFR after reddening correction. 
The submm brig htest phase has presumably a short d uration 
(~ 40 — 200Mvr: ISmail et al.l20()4tlGreve et alJ2005^ ■ while 
the subsequent normal sBzK phase is estimated to last quite 
longer (~ 0.5-1 Gyr; Daddi et al. 2005), which would explain 
why SMGs have a factor of 10 or more lower spatial density 
llScott et a l J2002tlDaddi et al.^2004;j bji. while havin g similar 
masses te.g.- lGenzel et al .l uOOSn^eve et alJl2005h a nd per¬ 
haps spatial clustering llBlain et alJl2004 iDaddi etani2004al: 
lAdelberger et al .1120051) . This picture is consistent with the 

evolutionary scenario proposed by Vega et al. (2005). Chap¬ 
man et al. (2005) also suggest that SMGs evolve into near-IR 
bright galaxies. 

We conclude by suggesting that the follow-up of the most 
UV active sBzK could be a promising alternative for find¬ 
ing submm/mm sources, possible c omplementary to the r adio 
pre-selection method presented by lChaoman et alJ J2001I) for 
a similar redshift range z~ 1.5 — 2.8, or the Spitzer-fMIPS 
pre-selection (Lutz et al 2005). 
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